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ABSTRACT 

Physical models of circumstellar disks surrounding pre-main sequence stars are currently constrained 
primarily by the spectral energy distribution (SED) of the system. Here we present first results from 
a survey of nearby Herbig Ae stars using nulling interferometry. We measure the spatial extent of the 
mid-IR emission for the stars HD 150193, HD 163296, and HD 179218. The results indicate the size of 
mid-infrared emission around Herbig Ae stars may be much smaller than current models predict. The 
observations constrain ~90% of the 10 /iia flux to be within an approximately 20 AU diameter region, 
providing significant spatial limits for disk models. These results demonstrate the scientific potential for 
nulling interferometry to study circumstellar disks with better spatial discernment than is possible with 
standard mid-infrared imaging. 



Subject headings: stars:pre-main sequence - 
techniques : interferomet ric 



stars:individual(HD150193, HD163296, HD 179218) 



> 

in 
in 

O 



^ ■ 
Q-f 

6 ■ 

a: 



1. INTRODUCTION 

for emulateapj 

Circumstellar dust disks, the presumed nursery for plan- 
etary systems, are most easily studied while they are still 
young and are best seen around the brightest pre-main se- 
quence stars. Herbig Ae/Be stars are known to harbor an 
excess of infrared (IR) emission above that of their photo- 
spheres, and recent evidence has pointed to this excess be- 
ing in the form of a circumstellar disk (Natta et al. 2001). 
However, significant differences have been highlighted be- 
tween the Ae and Be spectral types indicating that while 
Herbig Ae stars do indeed harbor disks, as evidenced by 
detection of their millimeter flux, Herbig Be stars (earlier 
than B8), in general, do not (Natta et al. 2000). Herbig Ae 
stars are also interesting as the evolutionary precursors of 
the Vega-excess stars first discovered by the IRAS satellite 
to have infrared excess due to a debris disk. 

The spectral energy distribution (SED) has been a pow- 
erful tool for probing these circumstellar disks. Simple, 
geometrically thin, optically thick models of a disk fail to 
reproduce the observed spectral energy distribution, un- 
derestimating the amount of flux in the mid and far-IR. 
Recent models (Kenyon and Hartmann 1987; Chiang and 
Goldreich 1997; DAlessio et al. 1998) that can reproduce 
the SED predict a significant flare to the radial structure 
of the disk due to the vertical hydrostatic equilibrium, in- 
creasing the disk's exposure to the starlight. Absorption of 
this stellar flux heats the surface layer of the flared disk to 
a higher temperature than the interior. In addition to cre- 
ating emission features, the mid-IR flux from such systems 
is expected to be significantly more extended, compared to 
geometrically thin disk models. 

In this paper we present spatial constraints on disks 
around three Herbig Ae stars using a new technique called 
nulling interferometry. 

1 Observations reported here were obtained at the MMT Observa- 
tory, a joint facility of the University of Arizona and the Smithsonian 
Institution 



2. NULLING INTERFEROMETRY 

Nulling interferometry is a novel technique for detecting 
faint circumstellar structure around bright point sources. 
The technique overlaps the pupils of two telescopes or sub- 
apertures of a single telescope to interferometrically sup- 
press the light from an on-axis star. Light from an angular 
distance of = \/(2b) away will constructively interfere, 
allowing detection of faint circumstellar structure. The 
transmission pattern of the interferometer is given by 

T W =«dn»(^) (l) 

where b is the baseline of the interferometer, A is the wave- 
length of observation, and is the angle, along the line 
connecting the apertures, from the pointing center of the 
interferometer. Any extended object is seen in the focal 
plane of the interferometer as the object multiplied by the 
transmission pattern in equation (1) and convolved with 
the point-spread function of a single aperture. 

Our nulling interferometer is designed for use with the 
converted 6.5 m MMT, based upon a prototype first 
demonstrated with the old Multiple Mirror Telescope (Hinz 
et al. 1998). The nulling instrument, called the Bracewell 
Infrared Nulling Cryostat (BLINC; Hinz et al. 2001), di- 
vides the pupil of the telescope into two halves and over- 
laps them on a 50% transmissive beam-splitter. The mid- 
infrared array camera (MIRAC; Hoffmann et al. 1998) 
provides a final stop for the overlapped beams defining 
two 2.5 m apertures separated by 4 m. At 10.3 fim this 
provides a constructive fringe 0.27" from the star. The 
orientation of the interferometer can be adjusted since the 
MMT is an altitude-azimuth mounted telescope equipped 
with the ability to derotate the instrument. The derotator 
stays fixed for nulling observations, with its orientation de- 
termining the interferometer baseline's position angle on 
the sky. For a more detailed description see Hinz et al. 
(2001)'. 

Precision nulling interferometry is designed for use with 
adaptive optics which corrects the atmospherically aber- 
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rated wavefront to improve the destructive interference be- 
tween the beams. However, useful observations can be 
taken for bright objects even without atmospheric cor- 
rection. When the interferometer is properly aligned and 
phased on an object the atmospheric phase variations be- 
tween the apertures causes the image to randomly fluc- 
tuate in brightness. In effect the atmosphere acts as a 
phase shifter for the interferometer. By recording short 
exposure (0.05-0.5 s) images of the object we accumulate 
frames which show the star both at constructive and de- 
structive interference. If a sequence of frames is recorded 
on the object, the nulled intensity measured by the instru- 
ment (instrumental null) is defined as 

N mst = £™ (2) 

Imax 

where I m m and I m ax correspond to the minimum and max- 
imum intensity of a recorded sequence. By comparing this 
nulled intensity to that of a calibrator star, any excess 
flux above that of a point source can be reliably detected 
down to the level of uncertainty in the measured null. The 
measurements can also be thought of as measuring the vis- 
ibility of the object, which is given by the relation 

-it 1 Ni ns t , . 

Vinst — ~, j 77- • {O) 

In interferometry the true source visibility, given an in- 
strumental visibility (Vi) and a calibration point source 
visibility (V c ), is defined as V s = y-. From this it can be 
shown that the true nulled flux, or source null, is given 
approximately as 

N s ~ N, - N c . (4) 

3. CIRCUMSTELLAR DISK MODELLING 

Models which have been able to successfully address the 
SED of T Tauri and Herbig Ae stars are dominated by 
reprocessing of the stellar radiation via a flared disk. This 
creates a surface layer which is superheated by the stel- 
lar irradiation (Kenyon and Hartmann 1987; Chiang and 
Goldreich 1997; D'Alessio ct al. 1998). For the purpose of 
estimating the size of mid-IR emission in such a disk we 
use the model of Chiang and Goldreich (1997, hereafter 
CG). 

The spatial distribution of the IR emission due to the 
heated atmosphere of the disk is determined by the tem- 
perature distribution as well as the optical depth of the 




Fig. 1. — Sketch of circumstellar disk model (Chiang and Gol- 
dreich 1997). The star irradiates the disk at grazing incidence pen- 
etrating to an optical depth of unity and creating a layer of dust 
which is heated above that of the disk interior. The optical depth 
of the heated surface is proportional to the incidence angle of the 
starlight. 



atmosphere. CG derive a temperature distribution for the 
dust atmosphere of 

(R \ 2/s 

T atm {r) * T„ (^J (5) 

where T» and R* are the effective temperature and radius 
of the central star. If the stellar radiation strikes the disk 
surface at an angle, a, as shown in Figure 1, then, from 
geometrical considerations, it can be seen that the optical 
depth of the superheated layer is proportional to a, since 
the stellar radiation will penetrate to unity optical depth 
at grazing incidence. For a disk in hydrostatic cqulibrium 
the grazing angle is given by CG as 



a= 7 



4/7 / r N 2/7 
J?* 



(6) 



for large distances from the star, where T c — G ^f^ la , 
is the mean molecular weight of the gas, k is the Stefan- 
Boltzmann constant, and G is the gravitational constant. 
The optical depth is of the atmosphere is then given as r = 
a em, where cir is the infrared emissivity of the material. 

The interior of the disk also contributes a significant 
portion of the mid-IR flux. CG derive a temperature dis- 
tribution for the interior of the disk of 

™-.(?r(Tf- m 

which is valid for the region where the disk interior is 
opaque to its own radiation, true for mid-IR emission cal- 
culations. 

For an optically thin disk atmosphere and optically thick 
interior with temperatures as listed above, the emission 
contributed by a given annulus, dr, of a face-on disk is 



I v {r) = 2tt (r(r) B v {T atm ) + B V {T % )) r dr 



(8) 



where B„ is the Planck function. For stars with parameters 
similar to those listed in Table 1 the 10.3 /jm emission from 
the disk atmosphere is predicted to come from an annulus 
centered at ~17 AU from the star. 

The response of the interferometer is found by multiply- 
ing equation (8), scaled appropriately for the distance to 
the star, by the transmission pattern of the interferometer 
(equation (1)) to get 

'it bx s 
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dxdy. (9) 



\ v is only nonzero for values of r > 0.6 AU, corresponding 
to the sublimation region of silicate dust around A stars. 
The outer limit, r D , is taken to be 500 AU. The model is 
not sensitive to the exact value of these limits since little 
flux is produced at the extremes. Since the object size is 
on the order of the fringe spacing of the interferometer the 
light detected from the disk at constructive interference 
must also be calculated. This is given by 

hon = 4 jf ° jf ° J„(r) cos 2 (^^j dxdy. (10) 

The calculated expected nulls for the three Herbig Ae/Be 
stars are listed in Table 1. 



Nulling Observations of Hcrbig Ae Stars 
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4. OBSERVATIONS 

Observations were taken during the commissioning run 
of BLINC (June 10-18, 2000) at the MMT, shortly after 
first light of the telescope. Two sequences of 500 frames 
were recorded for each object at orthogonal orientations 
of the baseline (except for HD 150 193 which was only ob- 
tained in one orientation). An integration time of 0.5 s 
was used for each frame. All objects were observed June 
14, 2000 during a night of particularly good seeing. Ob- 
servations were taken using a 10.3 /im, 10% bandwidth 
filter. The images are sky subtracted using frames taken 
of adjacent sky between observations. A custom IRAF 
aperture photometry program is used to extract the min- 
imum and maximum frames in a sequence to determine 
the instrumental null. Error in measuring the instrumen- 
tal null is estimated by comparing subsequent sequences 
on the same object. Since the object can be detected even 
in the nulled image the error in estimating the null is dom- 
inated by the atmospheric variations in phase. Using this 
approach, measurements of the instrumental null are typ- 
ically measured to a precision of several percent. 

Calibration of the observations was performed using a 
Herculi and Vega, using the same integration time as for 
the Herbig Ae stars. The interferometer was scanned in 
path-length (by moving the star across the array where an 
arcsecond corresponds to 38 ^m of path difference) after 
any significant telescope movement or change in orienta- 
tion which might change the relative paths of each arm 
and decrease the visibility. 

We also observed a Scorpii during the same night for 
a different science program. The star was measured to 
have a source null of 19% (visibility of 68%) at 10.3 /im 
due to the extended dust outflow surrounding the star. 
This measurment agrees with other mid-IR intcrferomctric 
measurements such as the Infrared Spatial Interferometer 
(Danchi et al. 1994) and demonstrates the instrument's 
ability to detect extended dust structure. 

5. DISCUSSION 

As shown in Table 2 the measured nulls for three nearby 
Herbig Ae stars are indistinguishable from the null mea- 
sured on Vega and a Herculi (our point source calibra- 
tions) under the same conditions. Although the level of 
nulling achievable is not dramatic, due to the presence 
of atmospheric turbulence, the measurements clearly do 



not agree with the model presented of mid-IR emission for 
Herbig Ae stars. 

Meeus et al. (2000) present ISO spectra for a sample of 
Hcrbig Ae/Be stars including the three objects observed 
with nulling. All three objects show significant silicate 
emission features over the 10.3 /im filter pass band used 
for MMT observations. The strength of emission is similar 
for the three stars and contributes approximately 50% of 
the total flux in each case. The relative equality of the two 
contributions requires that both sources of the flux, the 
silicate emission and continuum, be spatially unresolved, 
since either one, if extended, would contribute more flux 
to the null than is detected in the measurements. 

As pointed out by Bouwman et al. (2000), present mod- 
elling of circumstellar material to fit an object's SED is de- 
generate with respect to the size distribution of the grains 
as well as their geometric distribution. For example, larger 
grains closer in to the star can mimic the spectra of smaller 
grains further from the star, as long as the grain size distri- 
bution is consistent with the 10 /im emission feature seen 
by Meeus et al. (2001). A similar degeneracy exists for 
constraining the inner and outer boundaries of the dust 
region. Thus the observations presented are best viewed 
in terms of breaking the current degeneracies in models of 
circumstellar material using the SED alone, by providing 
spatial constraints for the mid-IR dust emission. 

One might think that the observations can be explained 
by nearly edge-on disks which shield our line of sight to 
the flared disk atmospheres. However the silicate emission 
seen by ISO suggests we are seeing optically thin irradiated 
dust, suggestive of a surface layer or envelope. In addition 
inclined models by Chiang and Goldreich (1999) predict 
a significantly different SED for inclinations greater than 
60°. For inclinations less than this, even the minor axis of 
the projected disk would be detectable with our interfer- 
ometer. 

Independent of interpretation, however, the observa- 
tions constrain the spatial extent of circumstellar disks 
around Herbig Ae stars. To fit the observations the con- 
volution of a model flux with the transmission pattern of 
the interferometer must be less than the listed values in 
Table 2. Roughly, this corresponds to 90% of the 10 /im 
flux of the disk coming from within 0.13" diameter or 20 
AU at a distance of 150 pc. Observations of a larger sam- 
ple of Herbig Ae stars are needed to verify this result. 



Table 1 



Parameters and Expected Source Nulls of Herbig Ae Stars 



Star 


Spectral 


distance a 


rj-i b 






Expected 




Type 


(pc) 


(K) 


(R©) 


(M G ) 


Source Null 


HD150193 


A2Ve 


150 


10000 


2.1 


2.6 


63% 


HD163296 


AOVe 


122 


10500 


2.8 


3.3 


72% 


HD179218 


B9Ve 


240 


10000 


2.5 


3.4 


54% 



a Perryman et al. (1997). 
b Meeus et al. (2001). 
c Hillenbrand et al. (1992). 
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Of the three stars, the most detailed information about 
the surrounding excess is known for HD 163296. Mannings 
and Sargeant (1997) observed it at 2.7 mm and found an 
elongated dust disk, at an inclination of 58° approximately 
110 AU in diameter with a position angle of 126°. Grady et 
al. (2000) confirm this structure in the visible with HST, 
where a disk illuminated by scattered light is seen with an 
inclination of 60° and a position angle of 140°. The posi- 
tion angles corresponding to the two nulling observations 
for this star are 10° and 94°. Although neither of these 
lies along the major axis of the disk, the two observations 
would still be expected to resolve the disk at the interme- 
diate angles of observation if the mid-IR emission followed 
the spatial structure of the CG model. At an inclination 
of 60°, if the inclination does not cause obscuring of the 
emitting region, the expected source null would still be 
60% if the disk major axis were oriented perpendicular to 
the baseline direction of the interferometer. 

6. CONCLUSIONS 

The observations of three Herbig Ae stars using nulling 
interferometry constrains disks around these objects to be 
less than approximately 20 AU diameter in the mid-IR. 
Although more objects will need to be observed to confirm 
this conclusion, these first observations, if indicative of the 
sample, provide a mid-IR spatial constraint for models of 
proto-planetary disks around Herbig Ae stars. 

These observations demonstrate that nulling interfer- 
ometry is a useful technique for probing pre-main se- 
quence circumstellar environents in the thermal IR. The 
constraints placed on the spatial extent of the disk, even 
though they were achieved using a single aperture tele- 
scope, are more stringent than is possible by using stan- 
dard imaging. Typical widths of images of the models 
presented in §3 would result in an increase of roughly 20% 
in the object's FWHM, compared to a point source. While 
nominally measurable, effects such as chopper jitter, atmo- 
spheric variations, and variable optical aberrations could 
all work to mimic such a small increase in image width. 
In contrast the interferometer converts the observation to 
a relative photometric measurement, which is straightfor- 
ward to measure and calibrate. 

We would like to thank Marc Kuchner for suggesting 
the attractiveness of using nulling to probe Herbig Ae/Be 



Table 2 

Results of Nulling Observations 



Star 


Instrumental 


Source a 


Position 




Null 


Null 


Angle 


HD150193 


13±5% 


0±5% 


97° 


HD163296 


12±7% 


-1±7% 


94° 


HD163296 


16±2% 


3±3% 


10° 


HD179218 


16±2% 


3±3% 


162° 


HD179218 


14±2% 


1±3% 


87° 



a Calibrator stars were a Her and Vega with measured nulls 
of 13±2% and 13±3% respectively using the same integration 
time as the science objects. 
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